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Dark matter pair production at high energy colliders may leave observable signatures in 
the energy and momentum spectra of the objects recoiling against the dark matter. We 
use LEP data on mono-photon events with large missing energy to constrain the coupling 
of dark matter to electrons. Within a large class of models, our limits are complementary 
to and competitive with limits on dark matter annihilation and on WIMP-nucleon scatter- 
ing from indirect and direct searches. Our limits, however, do not suffer from systematic 
and astrophysical uncertainties associated with direct and indirect limits. For example, we 
are able to rule out light (< 10 GeV) thermal relic dark matter with universal couplings 
exclusively to charged leptons. In addition, for dark matter mass below about 80 GeV, 
LEP limits are stronger than Fermi constraints on annihilation into charged leptons in dwarf 
spheroidal galaxies. Within its kinematic reach, LEP also provides the strongest constraints 
on the spin-dependent direct detection cross section in models with universal couplings to 
both quarks and leptons. In such models the strongest limit is also set on spin independent 
scattering for dark matter masses below ~ 4 GeV. Throughout our discussion, we consider 
both low energy effective theories of dark matter, as well as several motivated renormalizable 
scenarios involving light mediators. 
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1. INTRODUCTION 

The search for dark matter and its interaction with standard model particles is actively pursued 
by experiments worldwide. Direct detection searches look for a feeble kick that a dark matter 
particle produces in recoiling off a nucleus. Indirect searches aim at the detection of the annihilation 
products of dark matter particles with each other in regions with a high density of dark matter. A 
signal in any experiment using either of these techniques requires the existence of a new interaction 
between dark matter and standard model particles. Direct and indirect searches, together with 
assumptions on the astrophysical dark matter density and velocity distributions, place bounds on 
such possible interactions. 

The very same interactions may also lead to the production of dark matter at a high energy 
collider (with an appropriate beam of incoming particles). In this article we will explore possible 
couplings of dark matter to leptons and the limits on such couplings from the LEP experiments at 
CERN. There, the annihilation of an electron and a positron into an invisible dark matter pair may 
become visible if an additional hard photon is radiated during the collision, producing a distinct 
mono-photon signal. Since the LEP experiments did not observe an excess of mono-photon events 
beyond the expected background, a limit may be placed on the postulated interaction strength 
between dark matter and the standard model. These limits, in turn, can be reinterpreted as limits 
on both direct and indirect detection rates, independent of astrophysical and atomic uncertainties. 

Previous work relating collider searches to direct and indirect searches for dark matter has 
focused on the Tevatron |TJ [2] and the LHC [3] . While these hadronic machines probe the dark 
matter couplings to light quarks, the LEP data we are going to study is sensitive to the dark 
matter-electron coupling. The potential limits from ILC mono-photons on a thermal relic that 
couples to leptons was studied in [1] . If dark matter were hadrophobic, as has been discussed j5j- 
[7] (but disfavored [U [9]) as a possible explanation of the DAMA [10] and CoGeNT [H] signals, 
as well as various cosmic ray anomalies, the LEP mono-photon searches would provide the only 
sensitive, model independent, collider limits for dark matter. As we shall see, LEP searches can 
yield bounds on dark matter which are both competitive with, and complementary to, those placed 
by traditional dark matter searches. 

The plan of this paper is as follows. In the next section we will introduce the effective theory 
formalism we will use in the first part of the paper. The list of operators we are going to consider will 
not be exhaustive, but will encompass the phenomenologically most relevant scenarios. We include 
cases where the dark matter-lepton couplings are scalar, vector and axial-vector in nature 1 which 
covers a broad range of phenomena, including spin independent and spin dependent scattering as 
well as annihilations which are either velocity suppressed or not. In section [3] we will set limits 
on the various contact operators from the mono-photon search at LEP. Then, in sections [4] and [5] 
we will translate our limits into bounds on dark matter nucleon scattering and dark matter self 
annihilation, respectively. We will compare our results to current direct and indirect searches. 
In section [6] we will consider the possibility that the effective theory described in section [2] is not 
appropriate for calculating the production rate of dark matter pairs at LEP. We will discuss several 
renormalizable models in which a new gauge boson or a new scalar particle is introduced to mediate 
the interactions of dark matter with leptons. As we shall see, the inclusion of such particles can 
significantly alter LEP bounds, and in certain regimes the bounds become sensitive to the details 
of the UV completion. We will conclude in section [7j 



1 Throughout we consider the dark matter to be a Dirac fermion, since our bounds would not be altered significantly 
if dark matter is a Majorana fermion [THS]. We also do not consider scalar or vector dark matter, though we do 
not expect the limits to be qualitatively different. 
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2. THE INTERACTION OF DARK MATTER WITH LEPTONS 

In order to produce dark matter at LEP it must couple to electrons. In many models this may 
occur via the exchange of a heavy mediator that can be integrated out of the theory at low energies. 
In that case one can describe the phenomenology in an effective field theory with higher dimension 
operators coupling the dark matter particle x to standard model leptons I = e, fi, r. This allows 
us to consider a large variety of dark matter phenomena without committing to a particular high 
energy framework 2 . We will be considering the operators 

(vector, s-channel) (1) 

(scalar, s-channel) (2) 

(axial vector, s-channel) (3) 

(scalar, i-channel) (4) 

which capture the essential dark matter and collider phenomenology (e.g. spin dependent and spin 
independent scattering on nucleons as well as s- and p- wave annihilation). The classification of 
these operators as s-channel or t-channel refers to their possible UV-completion: are most 

straightforwardly obtained in models in which dark matter is produced at LEP through a neutral 
s-channel mediator, while eq. Q arises most naturally if the mediator is a charged scalar exchanged 
in the t-channel. With such a UV completion in mind, the suppression scale A can be interpreted 
as the mass of the mediator M, divided by the geometric mean of its couplings to leptons, gi, and 
dark matter, g x : A = M / \jgng x . Note that we assume lepton flavor to be conserved in the dark 
matter interaction. LEP can only constrain couplings to electrons, £ = e, and in principle the 
suppression scale A could be different for couplings to fi and r leptons. In the following discussion, 
we will therefore consider both scenarios in which dark matter couples only to electrons (i.e. A = oo 
for £ = /x, r) and scenarios in which dark matter couples in a flavor-universal way to all standard 
model leptons. Note that the last operator, eq. Q, may be transformed into a linear combination 
of the first three operators, plus pseudoscalar and tensor contributions, using the Fierz identities, 
but we include it separately here because it is a common outcome of supersymmetric theories. 

The effective theory described by equations Q-Q is always a valid description of processes 
with low momentum transfer, in particular dark matter-nucleon scattering in direct detection 
experiments. In high energy processes such as dark matter production at LEP or dark matter 
annihilation, the effective theory breaks down if the 4-momentum transfer is comparable to or 
larger than the mass of the particle mediating the interaction. In the first part of our analysis in 
sections [3}{5j we assume that this is not the case, and derive bounds on the operators Q-Q from 
LEP mono-photon searches, which we will then translate into constraints on direct and indirect 
dark matter detection cross sections. In section [6] we will investigate how these bounds change if 
the mediator of dark matter interactions is light so that an effective theory description is no longer 
possible. 
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2 Indeed, several recent studies have used effective theories to analyze and draw connections among dark matter 
experiments [12H16] . 
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Figure 1: Distribution of normalized photon energy in single-photon events at DELPHI. The agreement 
between the data (black dots with error bars) and both the full DELPHI Monte Carlo (solid yellow/light 
gray shaded histogram) as well as our CompHEP simulation (dotted histogram) is excellent. The blue 
shaded histogram shows what a hypothetical Dark Matter signal from e + e~ — > j\x would look like. We 
have assumed vector-type contact interactions between electrons and dark matter, m x — 10 GeV, and 
A = 300 GeV, see eq. ([!]). The peak at x 7 ~ 0.8 corresponds to the process e + e~ — > *yZ° — > jup, with an 
on-shell Z°. 



3. LEP LIMITS ON THE EFFECTIVE DARK MATTER ELECTRON COUPLING 

In this section we will consider the operators Q-Q and derive limits on their suppression scale 
A from mono-photon searches at LEP. While all four LEP-detectors have studied single photon 
events [T7], we will here focus on data from the DELPHI experiment |18| I19j . for which we were 
best able to simulate the detector response. The data was taken at center of mass energies between 
180 GeV and 209 GeV, but since in the analysis the events are characterized only by the relative 
photon energy x 7 = -E 7 /-E'beam> we can make the simplifying assumption that all data was taken at 
an energy of 100 GeV per beam. We have checked that the error introduced by this approximation 
is small. For our Monte Carlo simulations, we use CompHEP [20|, 12 lj . which allows us to include 
the effect of initial state radiation (ISR) which we find to be non-negligible. For example, we are 
only able to reproduce the height and width of the on-shell Z° peak in the x 7 distribution for the 
background process e + e~ — > juu (cf. Figure [TJ if ISR is included. 

To analyze the event samples generated in CompHEP, we use a modified version of MadAnaly- 
sis [22] , in which we have implemented the analysis cuts and efficiencies of the DELPHI analysis as 
well as energy smearing according to the resolution of the DELPHI electromagnetic calorimeters. 
In doing so, we closely follow ref. |18j . 

In DELPHI, central photons with a polar angle 9 (with respect to the beam axis) in the range 
45° < 9 < 135° are detected in the High Density Projection Chamber (HPC) with a threshold 
of x 7 > 0.06. We assume the trigger efficiency for photons in the HPC to increase linearly from 
52% at E 7 = 6 GeV to 77% at 30 GeV, and then to 84% at 100 GeV. The trigger efficiency is 
multiplied by the efficiency of the subsequent analysis, which we assume to increase linearly from 
41% at 6 GeV to 78% at 80 GeV and above. 

For photons with 12° < 9 < 32°, detected in the Forward Electromagnetic Calorimeter (FEMC), 
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the threshold is x 7 > 0.1. The trigger efficiency increases linearly from 93% at 10 GeV to 100% 
at 15 GeV and above, and the analysis efficiency is the product of a linear function, increasing 
from 57% at 10 GeV to 75% at 100 GeV, and a constant 89%, with the first factor coming from 
the analysis cuts, and the second one describing the loss of events due to noise and machine 
backgrounds. In addition, we impose an energy dependent angular cut 6 > (28 — 80x 7 )°. 

Very forward photons (3.8° < 9 < 8°) give a signal in the Small Angle Tile Calorimeter (STIC), 
whose threshold is x 7 > 0.3, and we assume the efficiency to be 48%, based on the (incomplete) 
information given in [THJ . We again impose an energy dependent angular cut 6 > (9.2 — 9x 7 )°. 

The above, calorimeter specific, efficiencies are augmented by an additional 90% efficiency factor, 
applied to all photons. We found it necessary to introduce this overall efficiency factor to gain 
agreement in normalization between our simulations and the results of DELPHI. 

The relative energy resolution, a E /E, is 0.043 © 0.32/ 'VE in the HPC, 0.03 © 0.12//E 
Q.ll/E) in the FEMC, and 0.0152 © 0.135/\/£ in the STIC, where E is in units of GeV. Here 
© means that the different contributions to the energy resolution function are statistically inde- 
pendent. For example, we simulate the effect of finite energy resolution in the HPC by shifting 
the energy of each HPC photon by an amount 0.043 E ■ r% + 0.32 \[E ■ r2, where r\ and r2 are 
independent Gaussian random numbers. Since we find that with purely Gaussian energy smearing 
we are unable to reproduce the broad tails of the on-shell Z° peak in the x 7 distribution (Figure [l]), 
we impose an additional Lorentzian energy smearing with a width of 0.052 E. This is motivated 
by a fit to the calorimeter response to monoenergetic electrons, obtained from ref. |23j . 

We have verified our modeling of the DELPHI detector by simulating the energy distribution 
of single photons in the Standard Model. As demonstrated in Figure [TJ the agreement with the 
data (black dots with error bars) and with the DELPHI Monte Carlo simulation (solid yellow/light 
gray histogram) is excellent. Only in the very last bin (x^ > 1), the observed number of events 
is ~ 4<j higher than the prediction by both Monte Carlo simulations, probably due to imperfect 
modeling of the detector resolution function. We therefore omit this bin in the following analysis. A 
straightforward x 2 analysis then yields x 2 /dof = 21.5/19 for our simulation, and x 2 /dof = 20.6/19 
for the DELPHI Monte Carlo. 

When setting limits on dark matter properties, we use our own simulation only for the signal 
contribution, while the predicted backgrounds are taken from the DELPHI Monte Carlo. The blue 
shaded histogram in Figure [l] shows what a typical dark matter signal would look like for the case 
of operator Oy, with a dark matter mass of 10 GeV and with A = 300 GeV. Since most of the 
signal events are in the low-x 7 region, where SM backgrounds are only moderate, and since the 
spectral shape of the signal is different from that of the background, we expect good sensitivity to 
the dark matter-electron coupling A -1 . 

Indeed, a x 2 analysis yields limits on the cutoff scale A of order 250-500 GeV for dark matter 
masses m x < 80 GeV (see Figure [2j). In this mass range, our limits on dark matter-electron 
coupling are slightly better than the limits on dark matter-quark couplings derived from Tevatron 
mono-jet events [H[2]. The Tevatron limits, however, do not yet include spectral information, and 
they extend to dark matter masses of several hundred GeV, while LEP is completely insensitive 
t° m x ~ 90 GeV for kinematic reasons. The normalized photon energy distribution is similar in 
shape for all the operators considered. This leads to similar limits on the operators from eqs. ([!])- 
Q at low dark matter mass. Only the limit on the strength of the operator {x^)(^x)i e Q- @> is 
somewhat weaker. Using the Fierz identities this operator may be converted to a sum of other 
operators involving a product of a dark matter bilinear and a lepton bilinear. There is destructive 
interference between these operators leading to a smaller production cross section for mono-photon 
events and thus a weaker bound on the cutoff scale for this operator. When the dark matter mass 
m x exceeds ~ 30 GeV, the limits on different operators scale differently with m x since at this point 
the dark matter particles are produced closer to threshold and the detailed dependence of the cross 
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Figure 2: DELPHI lower limits on the cutoff scale A of the dark matter effective theory for the four operators 
eqs. Q-Q as a function of the dark matter mass. The wiggles in the plot are due to limited Monte Carlo 
statistics. 



section on the final state velocities becomes important. 



4. LIMITS ON THE DARK MATTER NUCLEON SCATTERING CROSS SECTION 

The next step is to translate the limits on A into constraints on the dark matter-nucleon 
scattering cross sections probed in direct detection experiments. Since LEP can only probe dark 
matter-electron couplings, while direct detection experiments are most sensitive to dark matter- 
quark couplings, this translation cannot be done in a completely model-independent way. We thus 
consider two extreme possibilities, one in which the dark matter couples with equal strength to 
quarks as it does to leptons, and another in which dark matter couples only to leptons without 
coupling to quarks at tree level. Limits on other models, in which the ratio of lepton and quark 
couplings is different (e.g. coupling proportional to B — L), may be easily derived from these two 
cases, as we shall see below. 

In order to compute the dark matter scattering cross section off a nucleon, N = p,n, through 
one of the operators in Q-Q, we need knowledge of the nucleon matrix elements (N\0\N). We 
use the values of these matrix elements presented in pQ, with the exception of (N\qq\N) in which 
we follow [23] but use an updated [25] value of the pion-nucleon sigma term S^at = 55 MeV. 3 
As mentioned earlier Ot can be converted from a "i-channel" operator to a sum of "s-channel" 
operators by use of Fierz identities. Due to the relative size of the nucleon matrix elements it is 
sufficient to keep only the scalar s-channel contribution, which has a coefficient 1/4. Thus, for 
equal cutoff scale A, the direct detection rate expected from the operator Ot is the same as that 
expected from 

First we assume that the coupling of dark matter to all SM fermions, and in particular to all 
flavors of quarks, is identical to its couplings to electrons. In this case, the LEP bound on A can be 
immediately converted into an upper bound on the rate expected at direct detection experiments. 
We show these bounds in Figure [3] and we see that the limits on spin-independent WIMP-nucleon 



3 Note however that recent lattice determinations [26H29] of the strange quark content of the nucleus are considerably 
lower. The effect on our bounds, assuming equal coupling to all fermions, is small. 
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Figure 3: DELPHI upper limits (thick lines) on the cross section for dark matter- nucleon scattering compared 
to results from direct detection experiments (thin lines and shaded regions). The left-hand plot is for spin- 
independent scattering, as would come from operators Os, Oy 1 O t , and the right is for spin-dependent 
scattering through operator Oa- The spin-independent limits of CDMS and XENON-100 are taken from 
Refs. .30, and [31], respectively. The spin-dependent limits of DAMA, XENON-10, PICASSO, COUPP 
and SIMPLE are taken from Refs. [10], [32], [33], [34] and [35], respectively. The DAMA and CoGeNT- 
allowed regions are based on our own fit [36] to the data from Refs. [10] and [11]. Following 37 , we have 
conservatively assumed large systematic uncertainties on the DAMA quenching factors: = 0.3 ±0.1 for 
sodium and q\ = 0.09 ± 0.03 for iodine. All limits are computed at the 90% confidence level, while the 
DAMA and CoGeNT allowed regions are shown at the 90% and 3er confidence levels. 

scattering (left-hand plot) are competitive with direct detection results only for very light dark 
matter, m x < 4 GeV. The direct detection experiments become insensitive to such light masses 
due to their energy threshold, whereas there is no such low mass threshold at LEP. The high 
mass cutoff at LEP is reflected in the rapid deterioration of the upper bound at m x ~ 90 GeV. 
The LEP bound also applies directly to inelastic dark matter [38] . since the splitting between the 
two dark matter states of ~ 100 keV is inconsequential to the kinematics at LEP. However, such 
models typically require considerably larger dark matter-nucleon cross sections than elastic dark 
matter, since the splitting allows only the high velocity fraction of the dark matter to scatter. Our 
bounds derived from LEP rule out the very highest scattering cross sections in the parameter space 
consistent with DAMA [36], but still leave the bulk of the parameter space allowed. 

For spin-dependent scattering we expect the LEP bounds to be more competitive since there is 
little variation in the bound on A between the operators responsible for spin-independent scatter- 
ing (Oy and Os) and spin-dependent scattering (Oa), whereas constraints from direct detection 
experiments are much weaker than in the spin-independent case. The reason for this is that, unlike 
spin-independent dark matter-nucleus scattering, spin-dependent scattering is not enhanced by a 
factor A 2 , where A is the nuclear mass number. These considerations are reflected in the right-hand 
plot of Figure [3] where the LEP limits surpass direct detection constraints for m x < 80 GeV at 
which point the phase space for dark matter production at LEP again starts to shrink. 

If dark matter does not couple to quarks at tree level, but only to leptons (for simplicity we 
assume the coupling to fi and r is the same as that to e, our conclusions are not significantly altered 
even if the coupling were only to electrons), the power of the LEP limits improves dramatically. 
The reason is that in this case, dark matter-quark scattering to which direct detection experiments 
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Figure 4: Diagram for vector-type dark matter-proton scattering at the one-loop level. 



are sensitive is only induced at the loop-level [8]. 4 The cross section for loop-induced dark matter- 
proton scattering through the diagram shown in Figure [4] is 



fl-loop — 



18 2 7T 3 A 4 

t=e,fi,r 



■[J2 f(q 2 ,m e )} 2 , (5) 



where a is the electromagnetic fine structure constant, \i v = m p m x / (m p + m x ) is the dark matter- 
proton reduced mass, and the loop factor f(q 2 ,mi) is given by 



f(q 2 ,m e ) = \ 




5q 2 + 12m? + Q{q 2 + 2m?)Wl ^arcoth Wl ^ -3g 2 log [mf/A 



l renl 



(6) 



We take the renormalization scale A ren to be equal to A. Moreover, we make the approximation 
that all the dark matter is moving at the local escape velocity, which we take to be v x = 500 km/sec, 
and that the momentum exchanged in the scattering is maximal, i.e. the scattering angle is 180° 
in the center-of- momentum frame. This will overestimate the rate of observed recoils at a direct 
detection experiment and will lead to a conservative upper bound. With these assumptions the 
four- momentum exchanged between the dark matter and the target nucleus is q 2 = —Afi 2 v x , where 
H is the invariant mass of the dark matter particle and the target nucleus. 

The bounds on dark matter-nucleon cross sections quoted by direct detection experiments are 
derived from the actually measured dark matter-nucleus cross sections under the assumption that 
the dark matter couples equally to protons and neutrons and that the cross section is independent 
of q 2 . Here, however, it only couples to protons and there is a q 2 dependence in the loop factor 
f(q 2 ,mi). Thus, to enable a straight comparison, we rescale the quoted bounds on a p by A 2 /Z 2 x 
(Yle m i)/J2e f(q 2 / m e)) 2 i with q 2 = — 4/A^; and we take A ren = 500 GeV, the result is only 
very weakly sensitive to this choice. Note that (J5J) and ^ are only approximations in the effective 
theory formalism. The exact form of the loop factor depends on the embedding of the effective 
theory into a complete renormalizable model. 

In Figure [5] we show the LEP bounds on dark matter in the absence of tree- level couplings to 
quarks. Since loop-induced dark matter-nucleon scattering is forbidden for axial- vector interactions 



4 Dark matter-electron scattering is irrelevant in all direct detection experiments including DAMA [Sj and Co- 
GeNT [Sj. Even though DAMA and CoGeNT would not reject bulk electron recoils as background, kinematics 
dictates that the recoil energy can only be above the detection threshold if the electron enters the interaction with 
an initial state momentum > 10 MeV. The probability for this is very small due to the fast drop-off of the electron 
wave functions at high momentum [5] [H |9] . 
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Figure 5: DELPHI upper limits on the cross section for spin-independent dark matter-nucleon scattering 
for the case of dark matter with tree level couplings only to electrons, but loop level couplings also to 
quarks, compared to results from the direct detection experiments DAMA [TU], CoGeNT [TT], CDMS [50] , 
and XENON-100 [31]. The DAMA and CoGeNT allowed regions are based on our own fit [36] to the data 
from refs. [101 E]. We conservatively assume (?Na = 0.3 ±0.1 and gi = 0.09 ±0.03 for the DAMA quenching 
factors. All limits are computed at the 90% confidence level, while the DAMA and CoGeNT allowed regions 
are shown at the 90% and 3a confidence levels. 



and suppressed by two loops for s-channel scalar interactions jS], we consider only the vector- type 
operator Oy and the scalar i-channel operator Ot. As before, we apply the Fierz identity to Ot to 
decompose the operator into a linear combination of s-channel operators, of which we keep only the 
vector contribution. As is apparent from Figure [5j an explanation of the DAMA and/or CoGeNT 
signal by a dominantly leptophilic dark matter candidate which couples to nuclei only through 
loops is ruled out by LEP. 

Here we only considered two benchmark cases, where dark matter couples universally to SM 
fermions and when it couples only to leptons. Constraining a more general theory with a particular 
ratio of quark to lepton couplings, R q /u is straightforward. In this more general case nuclear recoil 
proceeds via both mechanisms, direct couplings to quarks and via a lepton loop. The limit on 
the former may be obtained by rescaling the bounds of Figure 3 by R^^ , whereas the limit on the 
latter may be taken directly from Figure [5} Generically one of these limits will dominate the other 
over the full dark matter mass range, and the less constraining bound should be taken. 



5. LIMITS ON THE DARK MATTER ANNIHILATION CROSS SECTION 

The LEP constraints on the suppression scale A of the effective dark matter couplings can 
also be converted to an upper bound on the annihilation cross section of dark matter into an 
electron-positron pair. They can then be compared to results from astrophysical probes of dark 
matter annihilation. Moreover, if dark matter is a thermal relic and if annihilation into electrons 
and positrons is the dominant annihilation channel, a lower bound on the dark matter abundance 
in the universe can be derived. If dark matter has also other annihilation modes, this bound is 
weakened by a factor l/BR(xx ~~ * e + e~). 
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Figure 6: LEP upper limits on the dark matter annihilation cross section (av), assuming that dark matter 
production at LEP and dark matter annihilation as probed by astrophysical and cosmological observations 
can be described by contact operators. In the upper left panel, we show limits on the process XX ~* e + e~ 
(the only one that can be constrained model-independently by LEP), while in the other panels we have 
made the assumption that dark matter couples equally to all charged leptons. For the average dark matter 
velocity (v 2 ) we have assumed the value at freeze-out in the top panels, while the bottom left panel is for 
the Draco dwarf galaxy which has very small (w 2 )- In the bottom right panel we compare the LEP limit on 
the u-independent interactions, Oy and Of, to limits from a variety of astrophysical observations [5§l441) . 



In order to translate the LEP constraints on the coupling strength A -1 into limits on dark matter 
annihilation, we need to calculate the annihilation cross sections corresponding to the operators in 
equations Q-Q. For annihilation into a single single lepton flavor of mass me, they read 



1 / TnP 

asVml = S^rA 4 " V 1 "4 (mX " m ' ) ^ ' (7) 
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'•..-) "73377 1/ 1 - — ! 21C2n4 + m?)+ x ^ r 2 cl J , (8) 

,,,, . ., 8mJ - 22m 2 mj + 17m 4 9 \ 
' M '"'- 1 = "n^TTl/ ' " — I 2 ''"< + m 2_ m | J > ( 9 ) 

mj ( 2 (m x + mi) 2 (8m 2 — 16m x mg + llmf) 



= 192^ \/ 1 " mj [ 24{m * + mt) + mf^ *W ) • (10) 



Here, we have made an expansion up to second order in the relative velocity v Te \ of the annihilating 
particles. While v TC \ <C 1 in all relevant astrophysical and cosmological environments, its exact 
value ranges from v ie \ ~ 0.1 at the time of dark matter decoupling in the early universe (if dark 
matter is a thermal relic) to values of order v Te \ <> 10 -4 (less than 100 km/s) in dwarf galaxies (see 
Appendix [A]) . This large spread of relative velocities can have a large effect on annihilation rates 
for certain operators. Notably, annihilation through an s-channel scalar operator, ([7]) is suppressed 
by v 2 el , and annihilation through an s-channel axial vector operator, ^ is suppressed by v 2 el or by 
m 2 /m x compared to the other modes. The production cross section at LEP is not suppressed in 
either of these cases, giving our bounds on the suppressed modes a substantial relative advantage 
compared to indirect searches. However, we will see that even in cases where the annihilation rate 
is unsuppressed the LEP bounds are interesting and competitive for light dark matter. 

In Figure [6j we consider both annihilation in the early universe and annihilation in the Draco 
dwarf galaxy 5 and compare to the cross section required for a thermal relic ((crv) ~ 3 x 10 -26 cm 3 /s) 
and to several astrophysical bounds. Our most model-independent bounds, those on annihilation 
into e + e~, are shown in the top left panel of Figure[6j where we take (^ei) = 0-24, corresponding to 
thermal freeze-out. We see that, if the dark matter only annihilates to electron-positron pairs, the 
thermal relic cross section is ruled out by LEP at 90% C.L. if m x < 20 GeV for vector interactions, 
and if m x < 50 GeV for scalar and axial vector interactions. Thus, in order for such dark matter 
to be a thermal relic it must have additional annihilation modes. 

Even though model-independently LEP can only constrain the dark matter coupling to electrons 
and hence the annihilation cross section for the process XX ~^ in many models of dark 

matter the annihilation rate into electrons is either equal or not very far from that into fi and 
r. For example, in models of supersymmetry the annihilation rate into charged leptons is set 
by the slepton masses, which in many cases differ by less than an order of magnitude. In other 
models, such as universal extra dimensions the annihilation rates to electrons, muons and taus are 
identical. 6 In order to present our results we pick the simple benchmark in which the operator 
strengths are universal among charged leptons. Constraints on other models may be derived from 
this benchmark by the appropriate rescaling. Limits on this benchmark scenario are presented in 
the upper right, lower left, and lower right panels of Figure[6j Due to the strong dependence of (av) 
on the charged lepton mass for axial vector interactions, the limit on the combined cross section 
for annihilation into all charged lepton species becomes significantly stronger below the r threshold 
in this case. In the lower left panel of Figure [6| we compare the LEP limits to constraints from 
Fermi-LAT observations of the Draco dwarf galaxy in gamma rays |39]. For m x < 80 GeV, LEP 
is superior to Fermi for all annihilation operators considered here, especially for scalar interactions, 



5 We chose the Draco dwarf galaxy because it is the dwarf galaxy for which Fermi-LAT obtains the strongest bounds 
on dark matter annihilation [39] . 

6 In models of universal extra dimensions, dark matter is usually a vector particle, a case we are not considering in 
this work. 

7 The Fermi-LAT collaboration presented their results as limits on the annihilation mode \X ~^ A £+ A t_ i assuming 
that this is the only annihilation channel. We have reinterpreted these limits, assuming that the branching ratio for 
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annihilation in our galaxy. Since the dark matter velocity distribution, especially at the galactic 
center is very uncertain, we include only the predictions for annihilation through the operators 



on annihilation through the axial vector operator Oa or the scalar operator Os will be between 
the corresponding constraints at freeze-out and those from the Draco dwarf galaxy and thus much 
stronger than the limits on vector interactions. Comparing the LEP constraint to limits from 
astrophysical observations, we find that the LEP limit is superior to Fermi results on gamma rays 
from dwarf galaxies |39| and on the high energy e + e~ spectrum [30]. We also find that the excess in 
gamma rays at the galactic center which has recently been argued |42j to plausibly arise from dark 
matter annihilations into r leptons is also strongly constrained, if the annihilation proceeds into 
electron-positron pairs at a similar rate. In fact, in |42| it was argued that an equal annihilation 
rate into electrons is favored because it may potentially provide an explanation of the WMAP 

Haze (mug. 

Constraints on dark matter properties from both indirect and direct observations are sensitive 
to the abundance, and velocity distribution, of dark matter both locally, at the center of the galaxy 
and in sub-halos. There are considerable uncertainties in all these quantities (32] that effect the 
exclusion curves, and preferred regions in Figure [6} We emphasize that the LEP constraints do 
not suffer from these astrophysical uncertainties. 



So far we have worked in a regime where the dark matter is the only particle of the dark sector 
accessible at colliders [32] and as a result all couplings of dark matter to the standard model are 
through higher dimension contact operators. However, since LEP is a high energy machine, there is 
a possibility that the particle that is mediating the interaction of dark matter with electrons is light 
enough to cause significant deviations from the mono-photon rates and spectra predicted by the 
effective theory. These deviations will be most pronounced when the mediator is produced on-shell 
and then decays to a dark matter pair, but as we shall see, order one deviations are possible even 
without on-shell production. We therefore also consider LEP bounds for several renormalizable 
"UV completions" of our effective theory. 

Possible renormalizable theories that couple dark matter with the standard model fall into two 
general categories, which we will dub "s-channel" and "i-channel" mediators. In the first case the 
mediator is a neutral boson which has coupling vertices to e + e~ and to dark matter pairs. In this 
case the mediator may be almost arbitrarily light if its couplings with matter are sufficiently feeble. 
Of the operators we consider here, s-channel mediators give rise to operators of the form of 0-([3]) 
at low energies. In the second case dark matter is produced at colliders via a t-channel diagram, 
exchanging a charged mediator. The canonical example is supersymmetry where neutralino dark 
matter may be produced at LEP by exchanging a scalar selectron. At low energies this gives rise 
to the operator Q. Since the mediator is charged in this case, its mass should exceed about 110 
GeV to evade direct LEP bounds. 

In cases where the momentum flowing through the mediator in collider environments is of order 
the mediator mass M, the momentum-dependence of the propagator has to be taken into account. 



the At + /i~ mode is 1/3 and that the 7-ray production is equal for all lepton flavors. In reality this will not be true, 
in particular there will be additional hard photon production for the r final state. A reanalysis of Fermi-LAT data 
including gamma rays from annihilation channels other than could improve the limits by an O(l) factor. 




suppressed by the small lepton masses. Limits 



6. CONSTRAINTS ON THEORIES WITH LIGHT MEDIATORS 
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Figure 7: DELPHI lower limits on the cutoff scale A = M / y/g e gx °f the dark matter effective theory. Dashed 
lines have been computed under the assumption that the effective theory is valid up to LEP energies, whereas 
the dotted and solid lines are for cases where the mediator mass M is so small that the effective theory 
breaks down. Once the mediator can be produced on-shell, its width F becomes relevant, as demonstrated 
by the shaded regions. r m i n is the minimum allowed width of the mediator, where g e ~ g x — M/A, and 
r min > lO" 4 GeV. 



In particular, the amplitude will be proportional to 

where q is the 4-momentum carried by the mediator, g e (g x ) is the coupling of the mediator to 
electrons (dark matter) and T is the total width of the mediator. In the case of an s-channel 
mediator q 2 = s — lyfsE^ is positive, while in the i-channel case q 2 is negative and depends on 
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the relative momentum between the two dark matter particles. In the previous sections, where 
the massive mediator could effectively be integrated out, the higher dimension operators were 
suppressed by a scale A. For the light mediator the LEP constraints become bounds on the 
geometric mean of the couplings g e and g x , but for ease of comparison we can still formally define 

A.^, (12) 

and quote the bounds in terms of that quantity. 

If the mediator and the dark matter are light enough to be produced on-shell at LEP, the 
bounds become sensitive to the width of the mediator T. V in turn depends on g e , g x (and possibly 
on the couplings to other particles). Here, we will treat T as a free parameter, but we note that, 
for any given value of ^Jg e g x (the quantity constrained by LEP), we can derive a lower limit r m i n 
on r by noting that 



The first term comes from decay into dark matter, the second one from decay into electrons, and 
'. . . ' stands for possible additional decay modes. For fixed ^g e g x the width is minimized if g e ~ g x , 
and if e + e~ and XX are the only allowed decay modes. If the latter assumption is true we can 



also place an upper bound on F by setting g e = Air and g x = M 2 /(g e A 2 m ) in (13), where An m is 
minimum value of A allowed by LEP. In what follows we will take the mediator's width to be a free 
parameter and will consider the effects of r m ; n < T < 1 GeV. For dark matter coupling through 
a t-channel mediator, no resonant enhancement is possible, so the value of T is irrelevant in this 
case. 

The limits on A = M/^/g e g x for various interactions are presented in Figure g From (ll|), 
we can understand the behavior of the dashed and dotted lines in this figure. Consider first the 
s-channel case for m x > M/2: There is no possibility of resonant production, so the mediator 
width is unimportant. Comparing the cross sections for dark matter production at LEP in the 
contact operator and light mediator cases, we obtain 



da 



dE 1 



M 4 da 



light mediator (q 2 — M 2 ) 2 dE. 



. (14) 

contact op. 



with q 2 = s — 2^Ey. Thus, for M slightly below -y/s, there is partial cancellation between the 
q 2 and M 2 terms in the denominator, leading to an enhanced cross section and an improvement 
in the limit on A compared to the contact operator case. For even smaller M, this cancellation is 
smaller and we expect the bound on A to scale with M. This is confirmed by Figure [7j 

On the other hand, if 2m x < M < y/s, the process e + e~ — > 7XX can proceed through an 
on-shell mediator, which leads to a peak in the monophoton spectrum reflecting the kinematics 
of a 2 — > 2 scattering process. The absence of a strong peak in the DELPHI data, apart from 
standard model Z production, places a strong constraint on this scenario. The constraint depends 
sensitively on the width of the mediator and scales as 

Aocl/r 1 / 4 . (15) 



This can be understood if we note that the resonant cross section for production of the mediator 
together with a single photon contains a factor 
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Figure 8: DELPHI lower limits on the cross section for dark matter-nucleon scattering for different dark 
matter interaction models. As in Figure [7J from which the limits are derived, dashed lines correspond to 
a contact operator interaction between dark matter and electrons at LEP, while the solid and dotted lines 
are for interactions mediated by light particles. In the background, we show the constraints from the direct 
detection experiments XENON-100, CDMS, DAMA, and CoGeNT (upper left, upper right and lower right 
panels) and from DAMA, PICASSO, XENON-10, COUPP and SIMPLE (lower left panel), see fig. [3] for 
details. 



where E Tes is the energy of the peak in the monophoton spectrum. Integrating (16) over the 
photon energy £L, we find that the total cross section for on-shell production of the mediator is 
proportional to 1/A 4 r (times factors that do not depend on A or F), which explains equation (15). 
We have also confirmed the scaling of the bound on A with r -1/4 numerically. 

Going back to Figure [7] and comparing the limits on A obtained for different types of operators — 
scalar, vector, and axial vector — we find that they are all comparable. The t-channel case is 
similar to the case of the s-channel away from resonance, except that the negative q 2 causes the 
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Figure 9: LEP upper limits on the dark matter annihilation cross section (av) for different assumptions 
on the mass of the particle that mediates dark matter production and annihilation. We show limits only 
for the annihilation channel \X e + e~ , which is the only one that can be probed model-independently at 
LEP. If dark matter has several annihilation channels, these limits can be straightforwardly (but in a model- 
dependent way) translated into limits on the total annihilation cross section, as done in the upper right and 
bottom panels of Figure [6j As in Figure [7J from which the limits are derived, dashed lines correspond to a 
contact operator interaction between dark matter and electrons at LEP, while the solid and dotted lines are 
for interactions mediated by light particles. 



denominator in equation (14) to be always larger than the numerator, meaning that the bound on 
A is always weaker in the light mediator case than for the contact operator. Furthermore, in the 
t-channel case there is obviously no on shell production of a mediator at low dark matter mass. 

Even though the effective theory is not appropriate to describe production at LEP, it is still 
a good description of dark matter-nucleus scattering in direct detection experiments, where the 
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exchanged momentum is very low. The procedure of translating our bound into direct detection 
limits is identical to that of Section |4} We present these bounds in Figure [8| making the assumption 
that dark matter has equal couplings to all standard model quarks and leptons. 

For non-resonant dark matter production (s-channel with 2m x > M or £-channel), the presence 
of the light mediator in general severely weakens the LEP bounds on the direct detection cross 
section. As discussed below equation (14), however, there is a window of mediator masses where 
the bounds are marginally improved compared to the contact operator case. If the mediator 
can be produced on-shell and is sufficiently narrow, the bounds on direct detection rates are 
strengthened considerably. In this case, the LEP constraints cover the (low mass) DAMA and 
CoGeNT- favored regions; for the vector operator a narrow resonance even impacts the DAMA 
region around m x ~ 50GeV. However, we should emphasize again that these conclusions can be 
evaded if the coupling of dark matter to electrons is much smaller than its coupling to quarks. 

Finally, we carry out a similar analysis to Section [5] and compute the annihilation rate in the 
early universe in the case of a light mediator. We consider only the case where the mediator couples 
exclusively to electrons and dark matter. Figure [9] shows that the LEP constraints on dark matter 
annihilation in the early universe change significantly if the mediator is light. The sharp peaks 
that occur at m x (l + (v 2 )/2) ~ M/2 are due to resonant annihilation of dark matter, and the 
dips observed just above the peaks are due to the fact that resonant annihilation and the on-set 
of resonant production at LEP occur at slightly different values of m x . As in Figures [7] and [8j the 
width of the mediator is of crucial importance for m x < M/2. 



7. CONCLUSIONS 

Very little is known about the dark sector of particle physics. It is usually assumed that dark 
matter couples, to varying degrees, to all fermions in the standard model, and strong constraints 
have been placed on its coupling to quarks by direct and indirect detection experiments and by 
the Tevatron. However, it is possible that dark matter has no coupling to quarks or at least that 
couplings to leptons are dominant. In such a scenario, dark matter may be efficiently produced in 
collisions of electrons and positrons at LEP. Irrespective of whether dark matter is leptophilic or 
not, LEP is an additional probe of its properties, and in this paper we have studied what LEP can 
say about the dark sector. Unlike dedicated dark matter searches in direct and indirect detection 
experiments, our LEP bounds do not suffer from astrophysical or atomic uncertainties. 

One mode in which dark matter may be searched for at LEP, with relatively little model 
dependence, is its pair production in association with a hard photon. The LEP experiments have 
searched for anomalous mono-photon events in their data sets, but have found no discrepancy 
from the prediction of the standard model. Unlike at hadronic machines, at LEP the kinematics 
of the event can be completely determined allowing the standard model backgrounds to be more 
easily distinguished from dark matter production. We used the mono-photon spectrum from the 
DELPHI experiment to place bounds upon the properties of dark matter that couples to electrons, 
see Figure [2} In the first part of the paper, we worked in an effective theory framework, in 
which dark matter interactions are described by four-fermion contact operators, and we derived 
constraints on the suppression scale, A, of these operators. 

We applied the LEP bounds on electron-dark matter coupling to constrain both the direct 
detection cross section and the annihilation rate of dark matter. We considered both the case 
where dark matter couples equally to all leptons and a scenario in which dark matter couples 
equally to all standard model fermions. Not surprisingly, for the "leptophilic" scenario, where 
LEP is probing tree-level interactions but direct detection proceeds through a loop process, LEP 
bounds are highly competitive. In fact, the bounds presented here rule out the DAMA favored 
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region, excluding leptophilic dark matter as an explanation of the DAMA modulated events or the 
CoGeNT excess, see Figure [5j 

In the case of equal couplings to quarks and charged leptons, the LEP bounds are complementary 
to direct detection bounds on spin-independent dark matter, see Figure [3j They are weaker than 
existing direct detection bounds for dark matter mass m x larger than ~ 4 GeV, but for light dark 
matter, m x < 4 GeV, they are significantly stronger. For spin-dependent interactions, where direct 
detection constraints are relatively weak, LEP outperforms all other experiments up to its kinematic 
limit, m x < 80 GeV. LEP bounds are slightly stronger than those derived in pQ from Tevatron 
mono-jet searches, but do not extend to as high masses, and they depend on the assumption that 
dark matter has universal couplings to quarks and leptons. 

We have also used LEP bounds to constrain dark matter annihilation rates, both in the early 
universe and in present-day galaxies. Below the LEP kinematic limit the LEP constraints are 
highly competitive. In particular, for m x < 80 GeV, they are stronger than those coming from 
Fermi-LAT observations of dwarf galaxies and of the galactic center, see Figure [6j They also 
provide a non-trivial constraint on a model invoked recently to explain a possible 7-ray excess at 
the galactic center [IT] . 

In the second part of the paper, we have repeated our analysis for the case where the interaction 
between dark matter and electrons cannot be treated as a contact operator. We have "UV com- 
pleted" the theory by introducing a particle that mediates dark matter-standard model interactions 
and have investigated LEP constraints as a function of the mediator mass and width. We find that, 
as long as dark matter cannot be produced through an on-shell mediator at LEP, our constraints 
are generally weaker than in the contact operator case (except for a narrow range of mediator 
masses close to the kinematic threshold of on-shell production). If the mediator mass M is below 
the LEP center of mass energy, but larger than 2m x , dark matter can be produced resonantly. In 
this case, the LEP constraint depends strongly on the width T of the mediator — a model-dependent 
quantity — but if T is small enough, the LEP constraint on the dark matter-electron coupling can 
be significantly stronger than for the contact operator case, see Figures [7j [8] [9| 

As the hunt for dark matter continues and we probe the dark sector on several fronts, both 
indirectly, directly and at the Tevatron and the LHC it is amusing to discover that there are 
non-trivial constraints still to be found in now completed experiments. It seems that dark matter 
requires us to be students of history as well as physics. 
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Appendix A: Derivation of the average dark matter velocity in a dwarf galaxy 

In this appendix, we discuss the derivation of the average dark matter velocity in the Draco 
dwarf galaxy. We assume the radial distribution of dark matter in Draco to follow a Navarro- 
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Frank-White (NFW) profile [17], 

„2 



p( r ) = Ps-r J ^2> ( A1 ) 

r (r + r s ) A 



with scale radius r s = 2.09 kpc and scale density p s = 0.98 GeV/cm 3 [39] 
We then use the Eddington formula 



1 f £ IT 1 d 2 p 1 dp 



(A3) 



to translate p(r) into the velocity distribution f(r,v) at radius r. Here, Vl/(r) = —Gj^drM(r)/r 2 
is (minus) the gravitational potential at radius r, which is determined by the enclosed mass M(r) = 
f r drAirr 2 p(r), and £(r,v) = \P(r) — \v 2 is (minus) the dark matter energy per unit mass. The 
dark matter density p is treated as a function of rather than r here, which is well-defined if \l/(r) 
is a monotonic function of r; the NFW profile has this property. The resulting velocity distribution 
f(r, v) satisfies the normalization condition 

p(r) = 47r Jdv v 2 f(r, v) . (A4) 

The annihilation rate of dark matter is proportional to p 2 crv, thus to obtain (f 2 ), the average dark 
matter velocity in the Draco dwarf galaxy quoted in Section [5j we compute 

i poo pi 

(v 2 ) = — driirr 2 p 2 (r) dv 47ru 4 /(r, v) (A5) 

N Jo Jo 

with the normalization constant N = J °° dr 4irr 2 p 2 (r) . We find (v 2 ^ ~ (34.7 km/s) 2 , which 
corresponds to (v 2 el ) ~ (69.3 km/s) 2 . 
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